Formation of the simplest stable negative molecular ion H3 in interstellar medium 
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We present the theory of radiative association of atoms and molecules, and we apply it to the 
(H2-H[~) van der Waals complex. We discuss the possibility for the H;^ ion to be formed in the 
interstellar medium in an environment with abundant ionized molecular hydrogen Hj. The obser- 
vation of H;^ would also be a probe for the presence of H~ in the interstellar medium. By computing 
the electronic structure of the H;^ ion, we determine its dipole moment, bound states, rotational 
constants, predissociated vibrational resonances and their lifetimes, and suggest a way to detect the 
ion in the interstellar medium. 

PACS numbers: 98.38.Dq, 95.30.Ft, 33.80.-b 



Many chemical reactions in the interstellar medium 
(ISM) are powered by cosmic rays: Atoms and molecules 
(mainly molecular hydrogen) are ionized by the radia- 
tion that provides sufficient energy to initiate a chain of 
chemical reactions in interstellar clouds leading to the 
synthesis of polyatomic molecules. A number of positive 
ions have been observed and identified in the ISM, in 
particular, the HJ^ ion. It is the simplest triatomic pos- 
itive ion that plays an important role in chemistry and 
evolution of interstellar clouds [3, Q, as its abundance 
is strongly related to the production of in the ISM. 
In contrast, only a few negative ions have been detected 
so far in the ISM: CsN", C4H-, CeH^, and CgH" [|. 
While quite stable, the simplest negative triatomic ion, 
H;^ (predicted to be bound by about 0.013 eV [1]) has 
not been detected so far in the ISM. In this article we ar- 
gue that the Hg ion is indeed formed in cold (below 150 
K) interstellar clouds, provided that (and free elec- 
trons) are available. Using a recently published potential 
surface Q, we determine its main spectroscopic proper- 
ties and suggest possible ways to observe it in the ISM. 
We analyze the possible formation mechanism of Hg" in 
the ISM in collisions between H~ and II2 . In addition as 
has only one bound state and, therefore cannot be 
directly observed, the detection of Hg would be a probe 
for the presence of H~ in the ISM, which is believed to 
exist in the ISM but has not been detected so far. 

The chemistry of interstellar clouds is initiated by ion- 
ization of molecular hydrogen by cosmic rays with a typ- 
ical rate constant C 3 x 10~^^s~^ in diffuse interstellar 
clouds [2]. (Cloud densities are ~ 10^ cm"'^ in diffuse 
and ~ 10^ cm~'^ in dense clouds.) The ionized molecu- 
lar hydrogen quickly forms Hj^ in collisions with II2, 
with a rate constant ~ 2 x 10~^cm'^/s The escaped 
electron has a large kinetic energy and undergos many 
elastic collisions with environmental H2 before its rether- 
malization. Possible inelastic e~-|-H2 collisions will lead 
to vibrational excitation of II2, and to dissociative at- 
tachment (DA) e~-|-Il2 H+H^ for collision energies 
above the threshold at 3.7 eV, with a cross section of 
about gda ^ 10^^^ cm^ [1, 01- Therefore, because the 
DA reaction rate per one H2 molecule is larger than the 
rate of electron production ^, Ve(JDAn{ii2) > Cj the 



ions are produced in the ISM with a (binary) rate con- 
stant mainly determined by the product ^n(Il2). Note 
that the same rate constant also determines the produc- 
tion rate of H;^ in the ISM 0. 

While the collisions between H~ and II2 molecules 
have been studied both theoretically [3, Q and exper- 
imentally [1, the structure of the H;^ ion has been 
rarely studied in the past Q , and a single indirect obser- 
vation has been reported in laboratory plasmas TOj . The 
H;^ ion is well represented as a van der Waals complex 
(II2 • • • ) 4] . The molecule has several rovibrational 
states, bound by about 20-100 cm^^ for the lowest ones 
(see Table m. According to the present study, there are 
also a number of predissociated resonances, which can 
be described as excited rovibrational states (j, Vd) of H2 
perturbed by H" , coupled to the dissociation continuum 
^2{j' ,v'j)+ll~ with energy of the dimer state E{j',v'^) 
lower than E{j,Vd)- We found that the widths of the 
broadest resonances in the low energy spectrum (below 
4000 cm~^) are in the range of 0.2 — 1.5 cm~^, which 
corresponds to lifetimes of 3.5 — 26 ps. 

To form a bound Hg" molecule in H2+H^ collisions in 
the ISM two mechanisms are possible: three-body recom- 
bination (TBR) or radiative association (RA): 
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The decay of Hg" in diffuse clouds is determined by col- 
lisions with H;^, other positive ions, and by interstellar 
radiation. The TBR rate constant k^b can roughly be 
estimated considering just geometrical (van der Waals) 
sizes of reactants and procedure outlined in Ref. pd| . 
We obtained fcsb - lO"^^ _ iQ-30 cmVs. With typi- 
cal number densities n = lO'^ — 10'* cm^^^ in the ISM, 
the three-body recombination as a way to form Hg" is 
expected to be much slower than possible two-body pro- 
cesses involving H~. Using the theory developed below, 
we estimated the RA rate coefficient /cra to be about 
4 X 10^^* cm^/s, corresponding to a survival lifetime of 
H^ ions tra = l/(fcRA"-) with respect to RA. The life- 
time should be compared with lifetimes with respect to 
other processes that remove H~ from the interstellar gas: 
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0,0,0,0,0,4i 


-105.01 


211.2 
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345. 
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-43.90 


345.1 
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-17.86 


187.5 
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0,0,0, 1,0, 


-17.7 


188. 
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0,1,0,1,0,^^ 


73.24 


311.9 


2.36 



TABLE L Computed binding energies (with respect to the 
lowest dissociation limit H2(0,0)+H~) and rotational con- 
stants of bound rovibrational states of H;^ relative to the Z 
and X axes. The quantum numbers are defined in the text. 



photodetachment of an electron from H~ and recombi- 
nation with positive ions. We have estimated the hfe- 
time with respect to photodetachment using the avail- 
able data about the photoabsorption spectrum of 
[T^ and intensity of radiation in typical diffuse inter- 
stellar clouds [1] and obtained the upper limit for the 
lifetime tpd ^ 10^^ s. We have also estimated the life- 
time r — \/{n+k'^) with respect to recombination of 
with positive molecular ions in the ISM, in particular, 
with Hjl" 10~^cm^'^, the corresponding rate con- 

stant fc+ - 10" cm /s is estimated using geometrical 
sizes of and H^). Therefore, we predict that the 
H;^ ion is produced in cold (especially dense) interstellar 
clouds consisted of molecular hydrogen. We also point 
out that if the assumption that an important fraction of 
all electrons produced by ionization of H2 would form 
is correct, the destruction of Hjj" in diffuse clouds would 
be then determined not only by the dissociative recom- 
bination but also the recombination with H~ . Such a 
hypothesis would solve the enigma of the anomalously 
large column density of Hj^ in diffuse clouds . 

We have performed the numerical calculation of bound 
states, their rotational constants, predissociated reso- 
nances, and cross-sections within the C omp lete Nuclear 
Permutation Inversion (CNPI) group [ijl D-^h of the 
molecule, using the Hj^ potential surface from Ref. 
The details of the numerical procedure are given in Refs. 
[fl [Tst . We solved numerically the three-dimensional 
Schrodinger equation for the molecule in hyperspheri- 
cal coordinates, separating hyperangles from the hyper- 
radius and using the slow variable discretization. Figure 
[1] yields an insight into the hyperspherical adiabatic 
curves calculated for J = 0. Each curve at large hyper- 
radii is correlated with a H2 (wd, j)+II^ dissociation limit. 
The lowest bound states and resonances can be charac- 
terized by the approximate quantum numbers j, il, vt, Vd 
defined below. 

Symmetry and approximate quantum numbers for H3" . 
The H;^ molecule consists of atoms with three identical 
nuclei, described within the CNPI group D^h- However, 
due to the van der Waals nature of the H;^ ion in its low 
energy states 01, the study of the lowest rovibrational 
states of H;^ can be performed within the Molecular Sym- 
metry group (or MS group [l3} C2V , which is a subgroup 
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FIG. 1: (Color online) hyperspherical adiabatic curves for 
J = Q. Only A\ and A2 vibrational symmetries are shown. 
The curve of the E symmetry are very similar to the Ai and 
A2 curves for energies below 4000 cm~^ (due to the barrier 
for the proton exchange [3|), and are not displayed here for 
clarity. Horizontal lines indicate positions of several bound 
states (reported in Table |T]) and resonances. (Here and below, 
a.u. means atomic units) 

of the CNPI group above). Therefore, in order to iden- 
tify allowed II3" states for ortho- and para-configurations 
of the total nuclear spin / (/ = 3/2 and 1/2, respec- 
tively) we derive below approximate wave functions and 
quantum numbers of the lowest rovibrational states. 

The Hg molecule has a definite total angular momen- 
tum J and its projection m on the z-axis of the space- 
fixed coordinate system (SCS). The orientation of the 
molecular coordinate system (MCS) with respect to the 
SCS is given by the Euler angles a, (3 and 7. The MCS 
Z-axis connects the center of the II2 dimer with the nu- 
clei of H", and its X-axis is in the molecular plane. In 
the present model, we assume that the H2 dimer is char- 
acterized by the vibrational quantum number Vd and its 
angular momentum j with projection on the Z axis. Its 
orientation with respect to Z is given by the azimuthal 
angle 9. Assuming that the vibrational and rotational 
motions are uncoupled, the H;^ wave function is repre- 
sented by (omitting a normalization factor for simplic- 
ity): 

\J,3,n,vuVd) = [Di^{a,l3,^)]* Pf{cose)\vd)\vt) , (3) 

where the associated Legendre polynomial Pj^{9) de- 
scribes the H2 rotational state, D:^Q{a, (3, 7) is a Wigner 
function, and Vt is a quantum of motion along the Z 
axis. The H;^ molecule is then considered as a symmet- 
ric rotor, with two of the three moments of inertia (along 
the X and Y axis) almost equal to each other, m does 
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not influence the energy of the state, and will be omit- 
ted. The energy of the state E{\J, j,n,Vt,Vci)) is given 
(approximately) by: 

Ei\J,j, vt,vd)) = + + BzjU + 1) + 

[BxJiJ+l) + iBz-Bx)n^] . (4) 

where Bx and Bz are the rotational constants with re- 
spect to the X- and Z-axes, respectively. The wave func- 
tions \J, j,fl,Vt,Vd) transform under C2v operators (12) 
(permutation of nuclei of the H2 dimer) and E* (inver- 
sion [l3|) as: 



{i2)\j,j,n,vt,vd) = 

E*\j,j,n,vt,va) - (- 



(5) 



where we assumed that H2 is in the X^J^'^ ground state 
and the trimer in the state. The final step is to en- 
sure that the approximate wave functions of have a 
proper symmetrization with respect to exchange of iden- 
tical nuclei, using the projection operators i-^ on a par- 
ticular irreducible representation (irrep.) F of the CNPI 
group D^h. The symmetrized states will be referred to 
as I J, j, fl, Vt, Vd, r). For certain combination of quantum 
numbers J, j and fi, some of the projections Pr are zero. 
It means that the corresponding irrep. F is not allowed 
for this set of quantum numbers [l6j . The general rules 
are derived from Eqs. ([5]): both (even and odd, labelled 
as prime and double prime) parities are allowed for non- 
zero n. If rj = 0, the parity is given by (—!)''• The E 
irrep. is allowed for any combination of J, j, Q (assuming 
that the parity E' or E" is given by the above rule). The 
A[ and A'( (resp. A'2 and A2) irreps. are allowed for even 
(resp. odd) j. 

With such wave functions, we can determine the al- 
lowed states for para- and ortho-H;^. Because the three 
nuclei are identical fermions, the total wave function (in- 
cluding the nuclear spin factor) of H;^ can only be of A2 
or A2 irrep. Since the para-Hg" nuclear spin part of the 
wave function transforms as E' in D^h, the allowed spa- 
tial irrep. could be E' or E" . For ortho-Hj^ {A[ irrep. of 
the nuclear spin), only A2 and A2 rovibrational wave 
functions are allowed. The lowest rovibrational state 
\0,0,0,0,0,A[) is not allowed for H;^ (but allowed for 
Dg"). The lowest allowed state is |0, 0, 0, 0, 0, £;') is the 
para-Hj^ state. The lowest A2 rotational state (lowest 
ortho-H;^) is 10,1,0,0,0,^^}. 

Theory of radiative association in dimer- atom colli- 
sions. In order to estimate the cross-section and the rate 
coefficient for RA, we develop a theoretical framework to 
treat the radiative association of a dimer and an atom. 
Our approach is based on theory developed by Herzberg 
[l3|, and later used by several authors [18, 19, 20] for di- 
atomic molecules, and for photoassociation of cold atoms 

In order to adapt the theory to triatomic systems, 
similarly to Ref. [lO], we express the Einstein coef- 
ficient Agi y'^.q E for the photon emission from an Hg 



rovibrational state specified by quantum numbers q = 
{J, j, ri, Vd, F} during a Il2-|-H^ collision at an energy E. 
After a photon of energy fko is emitted, the triatomic 
ion is in a state specified by the quantum numbers 
q' ~ {J' , j' ,v'd,T'} and by v'^ for the vibrational mo- 
tion of the II2 • ■ • H" van der Waals system. The Einstein 



coefficient Agiy^.g E is given (in a.u.) by: 



A„, 



q',v't;q,E — \rqiy^-q^E\ 



(6) 



where is the matrix element of the dipole mo- 

ment (with three components r'^, a = — 1,0, -|-1). The 
value \fq\v'i.;q,E\'^ cau be evaluated using a technique sim- 
ilar to the one presented in Ref. [l^: we used Eq. (14- 
33) of Ref. [1^ for the line strength and average it over 
the initial states. In the present model, each rotational 
state of H;^ is characterized by a single symmetric top 
rotational function, so Eq. ^ reduces to: 
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q',vi;q,E = ^i2J' + l)x 



A , , 

■^q'.Vf 

J2{j'M\<,Vd,^'\l^"\j,^,E,Vd,r) 



(7) 



J 1 J' 



where /i"^ is the a component of the dipole moment cal- 
culated in the MCS. The initial coUisional state in the 
above expression is energy normalized. The probabil- 
ity Pq' ,v'i.:q.E of an RA event is given by the Einstein 
coefficient divided with the current density in the flux 
of incident particles l/(27r) for the energy normalized 
wave function. Finally, the RA cross-section is given by 



T^Pq',v'^;q,E 



Ik' 



'^q',v'^:q,E 



Y,{j',^',<,v'd,T'\fi-\j,n,E,Vd,T) 



(8) 



J 1 J' 

n a -n' 



To obtain the cross section aq{E) for the formation of 
any H;^ bound state {q' , v'^}, we have to sum over q' and 
v'f. Since the nuclear spin is conserved during the RA 
process, there is no need to include the nuclear spin de- 
generacy factor: after averaging over the initial state and 
summing up over final states, the factor will be one. 

The rate constant is finally obtained by a standard in- 
tegration over a Maxwell-Boltzmann distribution. In the 
integration, the nuclear spin degeneracy factors {21 -1- 1) 
as well as the rovibrational energy of the initial state of 
H2 should be taken into account (see, for example, Ref. 
[2^ for details on the averaging procedure). The rate 
constant could be accurately calculated provided that the 
dipole moment functions are known over all the config- 
uration space, which is not the case with the available 
ab-initio calculations. Here we propose a rough estima- 
tion of the rate constant to verify if the RA process is 
competitive with other processes in the ISM leading to 
the removal of H~ (if present) from the interstellar gas. 
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For given values of f2 and f2', only one term in the sum 
of Eq. ([S]) is not zero. The 3j symbol and the symme- 
try = A'{ © A2 ® E" of the vector of dipole moment 
(/i^^, fjP, ^'^^) in the D^h group determine selection rules: 
J ^ J' = J±l; il — > SI' = f2, fizbl. In addition, parities 
of the initial and final states should be opposite. 

The largest vibrational dipole moment matrix element 
in Eq. ([5]) is expected when ^ — fl' (because |^°| is much 
larger than |At^^|), j = j' , and t^^ = v'^. When all internal 
"vibrational" quantum numbers {vd,^,j) are the same 
except Vt, the vibrational dipole moment matrix element 
can be estimated as /i*^ calculated at the equilibrium ge- 
ometry multiplied with the vibrational Franck-Condon 
overlap. The Franck-Condon overlap can be approxi- 
mated using the density of discrete states 1/AE calcu- 
lated at the energy of the final bound state v'^ . The value 
of 3j can be taken to be 1 for the allowed transitions 
( J ^ J ± 1; SI ^ 57' = ri) in the rough estimation. The 
equilibrium value of is 4 a.u. [l^], Ai? ^ 30 — 60 cm^^ 
or 1.5 — 3 X 10^** a.u.; fku ^ 100 cm^^ or 5 x 10"'* a.u.; 
k'^/{2m) ~ 1.3 X 10"'* a.u (it corresponds to 40K, a rea- 
sonable temperature for cold diffuse clouds), m ~ 1200 
a.u. is the reduced mass of the H2-|-H~ system. We took 
J = and, correspondingly, J' = 1. Plugging these val- 
ues into Eq. ([5]), we obtain the value of 10^^ a.u. for 
the estimated RA cross-section. The rate coefficient is 
estimated as cjq ■ k/m ~ 6 x 10^^^ a.u. or 4 x 10^^^ 
cm'^/s. In the estimation, we have neglected all the Fes- 
hbach resonances present in the coUisional spectrum of 
and H2. Such resonances should increase the total 



cross-section and will be accounted in an accurate calcu- 
lation in a separate publication. 

As a summary, the main results of the article are the 
following. (1) We have proposed a realistic scheme for 
the formation of the H3" ion in cold interstellar clouds. 
According to our results (2) the ion can be detected in the 
absorption spectrum in the millimeter wavelength range 
and, therefore, (3) can serve as a probe for the presence 
of in the ISM. To estimate the rate of formation 
in H2-l-H~ collisions, (4) we have developed a theoretical 
approach to calculate the cross-section for the radiative 
association and determined exact and approximate quan- 
tum numbers that can be used to characterize the bound 
and resonant states of H;^ . The developed theory is gen- 
eral and can be used to study other collision processes 
between a dimer and an atom, for instance, photoassocia- 
tion in ultra-cold gases. More details about the structure 
of Kg" and its isotopologues, including bound and reso- 
nant states will be presented in a separate publication, 
which will be relevant for considering new laboratory ex- 
periments with , for example, in ion traps 
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